Fermi-surface mapping from Compton profiles: Application to beryllium 
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The two-dimensional momentum density of Be on the basal TMK plane, i.e. the line integral 
of the three-dimensional momentum density along the c-axis, is reconstructed via the Cormack 
method from both experimental and theoretical Compton profiles. It is shown that in Be, despite 
the momentum density is highly anisotropic, merely two Compton profiles are sufficient to reproduce 
the main features of the momentum density. The analysis of the reconstructed densities is performed 
both in the extended and reduced zone schemes. 
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I. INTRODUCTION 

During the last few decades, Compton scattering has 
been widely used as a probe of the electronic structure 
of matter. Within certain approximations 1 —, the experi- 
mental double-differential inelastic x-ray scattering cross 
section can be related to the Compton profile, which is 
a double (plane) integral of the electron density in the 
extended momentum space, p(p), 



J(pz) = 



p(p) dp x dp y . 



(1) 



The Compton profile, being related directly to occupied 
electronic ground states, contains information not only 
on the Fermi surface (FS) but also on the Umklapp com- 
ponents of the electron wave functions. 

Eq. represents the Radon transform* of p(p) in 
three-dimensional (3D) space over 2D hyperplanes. The 
first reconstruction method based on the inversion of the 
Radon transform was Mijnarends' method 4 . Other alter- 
native solutions are Fourier transform algorithms^SiL^, 
the Maximum Entropy technique^, and two different 
methods employing orthogonal polynomialsiSiii. 

In cases where the electron momentum density is a 
strongly anisotropic function, reconstruction of the 3D 
momentum density from ID Compton profiles is rather 
elaborate, since a large number of measurements is 
needed. Thus a reconstruction of 2D momentum den- 
sities from Compton profilesi2ii£ is proposed, i.e. a con- 
version from double (plane) integrals in Eq. JIJ to a single 
(line) integral of the form 



much more easily than the Compton profiles. Further- 
more, it can be directly compared to corresponding 2D- 
ACAR (two-dimensional angular correlation of positron 
annihilation radiation) results, yielding additional infor- 
mation on the electron-positron interaction. It is also 
much easier to reconstruct accurately than the 3D mo- 
mentum density p(p) because the number of required 
Compton profiles is relatively small. 2D momentum den- 
sities can also be folded into the first Brillouin zone (BZ) 
in order to delineate various FS elements. 

This procedure is demonstrated for beryllium metal by 
using two Compton profiles, taken with the p z -axis along 
the reciprocal lattice vectors [100] and [110]. In this case 
the 2D momentum density can be reconstructed with the 
integration axis along [001]. The resulting densities are 
shown both in the extended and reduced zone schemes. 
For obtaining the density in the reduced zone scheme, 
the Lock-Crisp- West (LCW) procedure 14 is used. 



II. CONVERSION FROM ID TO 2D SPECTRA 

The purpose of this work is to reconstruct p^{j> z ,Py)i 
with L 1 1 [001], from its line integrals, i.e. from the Comp- 
ton profiles 



J(Pz) 



P^{Pz,Py)&Py 



(3) 



P^iPz.Py) = J(Pz,Py) 



P(p) &Px, 



(2) 



For this, the Cormack method— is proposed, corre- 
sponding to the solution of the Radon transform in 2D 
space over ID hyperplanes. 

If the direction p x is along an i?-fold rotation axis of the 
lattice, both the Compton profiles and the 2D momentum 
density can be expanded into a cosine series 



where L|[p x . When studying the structure of the Fermi 
surface, this 2D momentum density can be interpreted 



P^iPziPy) = P(p,@) 



]pf R (p)cos(tR&) (4) 



2 = 



J(Pz) = J{t,tp) = ^2gi R (t)cos(iR<p). 



(5) 



i=0 



Here (t,ip), where t = \p z \, are the polar coordinates of 
p z in the laboratory coordinate system, while (p, 9) arc 
the polar coordinates of the momentum (p z ,p y ) in the 
crystal coordinate system. 

By applying the Cormack method^, p L (p z ,p y ) can be 
evaluated from Eq. J3J| where its radial components p\ (p) 
for I = iR are given by 



Pi 



(6) 



fc=0 



R^{p) are the Zernike polynomials and the coefficients 
aik are obtained from the series expansion of gi(t) 




T 

■ First and H Second zone holes 
I 1 Third and g§§| Fourth zone electrons 



fji 



(t) = 2j2aik\/l-t 2 U l+2k (t), 



(7) 



k=0 



where U m (t) are the Chebyshev polynomials of the sec- 
ond kind. Since they are orthogonal in [—1, 1], 



1 



aik 



gi{t)Ui +2 k{t)dt. 



(8) 



To use the symmetry properties most profitably, the 
best choice is to measure J{p z ) for p z perpendicular 
to the [001] direction. The directions of p z (Eq. ©) 
should be equally spaced in the nonequivalent part of 
the BZ, which ranges in the TMK plane from ip = 0° 
up to 45° (or 30°) for cubic and tetragonal (or hexag- 
onal) structures, respectively. If N projections are to 
be measured for the structure with an i?-fold rotation 
axis along the [001] direction, the best orientation of 
the n-th projection (n=l,2,...,N) is given by the for- 
mula (p n = Aip/2 + (n — l)A<p where Aip = ir/(RN). 
However, if for any particular reason one would like to 
perform measurements also for the main symmetry di- 
rections, the orientation of the n-th projection ought to 
be determined from the formula ip n = (n — I) Aip' where 
Ay/ = n/[R(N - 1)]. 

In this way one can estimate J(p z ,p y ), which repre- 
sents the line integral of the density along the [001] di- 
rection. This procedure has already been applied to ID 
ACAR spectra of yttrium^. Even though the momen- 
tum density of Y is highly anisotropic, only three pro- 
jections were needed to obtain almost all details of the 
2D momentum density. In this work, the method is suc- 
cesfully applied to beryllium (another hep metal) in the 
extreme case of merely two measured projections avail- 
able. 



III. APPLICATION TO BERYLLIUM DATA 

The Be Compton profiles were measured at the Eu- 
ropean Synchrotron Radiation Facility (ESRF) beamline 



FIG. 1: Free-electron Fermi surface of Be in the TAK and 
TAM planes in the extended (upper panel) and reduced zone 
scheme (lower panel). 



ID15B. The details of the experiment are described in 
Ref. B Two projections were measured with p z along 
TM and TK and with an overall momentum resolution 
of 0.1 a.u. The corresponding highly accurate theoretical 
profiles were calculated with the KKR-LDA method 19 . 
Such spectra allow the reconstruction of p L = J(p z ,p v ) 
with p x along the hexagonal c-axis. This 2D momentum 
density will be denoted as p 001 (p z ,p y ). 

Let us choose the polar system (p' y ,p' z ), fixed to the 
lattice, with p' z along the TK direction. The measured 
profiles, J{p z ) = J{t, ip), are described by the polar angle 
ip = 0° for p z along TK and <p = 30° for p z along TM. In 
this case two radial functions gi(t) with 1=0 and 6 can 
be determined from Eq. (JSJ, 



go(t) = [ J rK (t) + JrM(t)}/2, 
56(0= I Jr K {t)- JvM{t)]/2. 



(9) 



It is convenient to choose the unit system so that t' = 
t/t max = cos(a;). The zeros of the Chebyshev poly- 
nomials U m (t') occur then at on = -£—(2% + 1) with 
i = 0, 1, . . . , m — 1, and 



M-l 



(l Ik 



= -^-{ V 3z(cos(«Aa)) sin[(2A; + I + T)iAa] 



1=1 



+ l(-V k+t/2 9i(0)}, (10) 

where Aa = n/2M , and M is the number of points used 
in the evaluation of aik. I n this work, M = 720 (Aa = 
0.125°), which allows one to calculate 150 coefficients aik 
for each gi (t). This is equivalent to applying the Gaussian 
quadratures for the polynomials up to the 300th order. 

The FS of Be for the free-electron model is presented 
in Figs. 1 and 2. 
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FIG. 2: The Fermi surface of Be in the TMK plane. Up- 
per panel: Free-electron model in the extended zone scheme. 
Lower left panel: Free-electron model in the reduced zone 
scheme. Lower right panel: Real FS derived from [14]. 
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FIG. 3: Differences between 2D densities of Be for momenta 
along TK and FM, reconstructed from two Compton profiles. 
The experimental results are indicated by solid squares. The 
corresponding results of the reconstruction from theoretical 
Compton profiles with and without resolution broadening are 
shown by open circles and solid line, respectively. 



To the authors' knowledge, the first realistic band- 
structure calculations for Be were performed with the 
OPW method 2 ^. The results of that study, compared to 
the free-electron model, suggested the following FS fea- 
tures: a) no holes around the H point either in the 1st or 
2nd band (1st zone fully occupied and no holes in the 2nd 
zone on the plane AHL); b) very small holes around £ 
and reduced holes around T in the 2nd zone when com- 
pared to the free-electron model; c) no electrons around 
r in the 3rd band; d) no electrons around L either in the 
3rd or 4th bands; and e) cigars in the 3rd zone around 
K are larger than for the free-electron model with their 
height close to \KH\ (see Fig. 2). These results were 
supported by de Haas-van Alphen (dHvA) experiment 21 
with two main differences: the computational cigars were 
triangular and the coronet was slightly larger than seen 
by the dHvA experiment. 

Further theoretical band structure results^Si were 
similar to those in [14] while the latest LCAO results^ 
turned out to be somewhat different, i.e. Fig. 3 in Rcf. 24 
representing p(p) along [Oil] = TL points out that there 
are electrons around the L point in the 3rd and 4th bands. 
Theoretical Compton spectra used in this workiS agree 
qualitatively with the band structure results in Ref. 20. 

The results in the extended zone scheme along two 
main symmetry directions TK and TM are displayed in 
Figs. 3 and 4. Fig. 3 presents the anisotropy of the 2D 
momentum density p 001 (p) for p || TK and p || TM, re- 
constructed from both the experimental and theoretical 
Compton profiles taken from Refs. llSlllPl The widely- 
used approximation to include the electron-electron (e-e) 
correlation into the momentum densities is the so-called 
Lam-Platzman (LP) correction— However, this correc- 
tion is isotropic and thus does not have any effect on the 
anisotropy of theoretical curves in Fig. 3. The resolution 



function of the experimental apparatus smears some of 
the fine structure of the data as indicated in the figure. 

For momenta 0.60 < p < 1.05 a.u. (atomic units of 
momentum), the reconstructed 2D densities have much 
larger values for p || TK than for p || TM. This is 
due to the fact that for the real FS there are no elec- 
trons in the 3rd or 4th bands around the L point. This 
is seen in Fig. [3] as a positive anisotropy. The largest 
positive values for \TM\ < p < \TK\ a.u. are, addition- 
ally, a manifestation of the lack of holes around the H 
point. The negative anisotropy for p > 1.05 a.u., with 
the minimum at p = p F M , originates from the fact that 
p T F M > pp , i.e. the holes around the £ point, shown 
in the lower panel in Fig. 2, are very small. Umklapp 
components for p > 1.9 a.u. along TK are also observed. 
These results are in good agreement with the theoretical 
3D momentum density^ and the experimental observa- 
tion of the Umklapp components 26 . All subtle features 
of the theoretical anisotropy are reproduced by the ex- 
periment in detail. The accuracy of this result is due 
to the fact that even with three Compton profiles avail- 
able (the third one measured with ip — 15°), the func- 
tion g 6 (t) and thus also Pq 01 (p) as well as the anisotropy, 
Pfk(p) — Pfm(p) ~ 2p6 01 (p), would be the same as ob- 
tained from two projections^. The fact that the mag- 
nitude of the anisotropy is diminished indicates that e-e 
correlation has an anisotropic effect on the momentum 
density, as has been suggested earlier—. 

Reconstructed densities for momenta along TK and 
TM are presented in Fig. 4. It can be seen that the ab- 
solute densities are not reproduced exactly, e.g. a small 
electron-like lens is observed at p — a.u. According 
to all previous band-structure calculations, this clement 
should not be present and is thus most probably an arti- 
fact originating from the fact that the isotropic compo- 
nent Pq 01 (p) is reconstructed from only two projections. 
However, this artifact is cancelled out in the difference 
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FIG. 4: 2D momentum densities of Be for momenta along TK 
(circles) and YM (triangles), reconstructed from two Comp- 
ton profiles. Theory (including resolution-broadening and 
the Lam-Platzman correction) and experiment are marked 
by solid and open symbols, respectively. 



shown in Fig. |31 which further demonstrates the accu- 
racy of the obtained anisotropy. 

In the analysis in the p space one should bear in mind 
that p 001 (Pz,Py) cannot be directly identified with the 
FS line dimensions along the [001] direction because in 
real metals p(p) < 1 also in the central FS. However, 
undoubtedly the theoretical as well as the experimental 
results show the lack of electrons around the L point 
in the 3rd and 4th bands and the shape of the FS on 
the TMA plane close to the double BZ boundaries. In 
order to obtain more detailed information, p 001 (j>z,Py) 
has to be folded to the reduced zone scheme (k space) as 
discussed in the following. 



IV. RESULTS OF THE LCW FOLDING 
PROCEDURE 

In the Compton-scattering experiment one measures 
the electron density in the p space 



occ poo 

p(p = k + G) = ]r| / ^{v)e-^dv\ 

J — OO 



(11) 



where k, G, j and 0kj(r) denote the wave vector, recipro- 
cal lattice vector, band index and electron wave function 
in the state kj, respectively. The summation is over all 
occupied states. 

In a periodic lattice potential, the momentum density 

in the j^ 1 band can be written as 



Pj (k + G)=n(kj) \u kj (G)\ 2 , 



(12) 



where u^j are the Fourier coefficients of the electron wave 
functions and n(kj) is the occupation number. 



It is well known that it is not possible to obtain the 
shape of the FS knowing only the density p(p). This is 
due to the fact that the density is not constant on the 
FS and p(p) represents a sum of contributions from all 
occupied bands, not only those crossing the Fermi energy. 
Thus, in order to obtain the FS map, the best choice is 
to perform LCW folding 14 , which is a conversion from p 
to k space. By doing this, the folded densities visualize 
the shape of the individual FS elements more clearly. 

In the case of the electron density (neglecting correla- 
tion effects), p|(k) = X^G_Pj(k+G) = 1 (° r for unoccu- 
pied states) and p e (k) = £ \ p|(k) = n, where n denotes 
the number of occupied bands contributing to the point 
k. The application of the LCW folding to p L (p z , p y ) gives 
the function p l '(k z , k y ) which represents the correspond- 
ing line integral of p(k), in this case along L = [001]. The 
density p L {k Zl k y ) can be identified with the sum of the 
FS line dimensions over occupied bands along L. 

The results of the folding procedure are displayed in 
Fig. 5, where densities p 001 (k z ,k y ) reconstructed from 
two theoretical and two experimental Compton profiles, 
are presented. The free-electron model result is also 
shown for comparison. 

Conventionally, the momentum densities p(k) obtained 
via the LCW folding, and thus also densities p 001 (k z , k y ), 
are presented in arbitrary scale. As was shown by Kaiser 
et al. |2^, the main reason is that the convergence of the 
LCW procedure is dependent on the degree of the local- 
ization of the electrons. The more localized they are, the 
larger the contribution of the Umklapp components at 
high momenta is, which consequently increases the trun- 
cation error of the LCW procedure. This effect decreases 
the LCW-folded densities in a nonuniform way, and thus 
an absolute normalization is difficult. In the case of ex- 
perimental densities, the finite resolution function and 
statistical noise, both affecting especially the Umklapp 
components, complicate the normalization even further. 

In Be the volume of occupied valence states is equal to 
the volume of two BZ's. Thus, if two bands were com- 
pletely filled, p 001 (fc 2 , k y ) would be constant and equal to 
2 |r^4| . The light and dark areas in Fig. 5 depict low and 
high electron densities originated from the holes in the 
second band (coronet) and the electrons around K in the 
third band (cigars), respectively. 

Fig. 6 shows the same densities in more detail along 
the main symmetry directions together with densities ob- 
tained from the corresponding LP-corrected theoretical 
Compton profiles. From Figs. 5 and 6 the 2nd-zone holes 
around T and the 3rd-zone electrons around K can be 
clearly recognized. Other features, e.g. the electron-like 
artifact around the T point in the 3rd band (seen also 
in Fig. 4) and a hole-like artifact between the S and M 
points are caused by using only two components of recon- 
structed density p 001 (p z ,p y ), determined from only two 
Compton profiles. 
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FIG. 5: 2D LOW densities p 001 (fc z , &„) in Be, in the repeated 
zone scheme, reconstructed from two ID projections. Panel 
(a): theory without LP correction, panel (b): theory with 
LP correction and resolution broadening, panel (c): experi- 
ment, and panel (d): free-electron model. Columns I and II 
show results described by 60 and 90 (or 120 in panel (a)) or- 
thogonal polynomials, respectively. Each figure contains 10 
contourlincs. 
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FIG. 6: 2D LCW densities in Be for momenta along TK 
(right-hand side) and TM (left-hand side), reconstructed 
from two ID spectra. The theory with LP correction with 
and without resolution broadening is indicated by solid and 
dashed lines, respectively, while the experiment is marked by 
circles. 



V. SUMMARY 

When the electron momentum density p(p) is highly 
anisotropic, it is quite difficult to obtain its shape from 
ID projections accurately, since a relatively large number 
of projections is needed. For these cases, the reconstruc- 
tion of 2D momentum densities J(p z ,Py) instead of 3D 
densities p(p) is advised, ft is shown that only a small 
number (2-5) of measured plane projections is needed to 
reconstruct the 2D momentum density, and generally this 
2D density allows much more detailed analysis of the FS 
than ID projections. The use of the Cormack method is 
proposed, because a) the expansion of measured spec- 
tra into orthogonal polynomials has mean-squares ap- 
proximation properties and thus it effectively smoothes 
the statistical noise in the experimental data2£; b) the 
method employs the Chebyshev polynomials, which are 
the only orthogonal polynomials having analytical zeros, 
allowing the estimation of the coefficients aik to a high 
degree of precision. 

This method was applied to two Be Compton profiles 
and the line projection of the momentum density along 
the [001] direction was reconstructed. 

The analysis performed in the extended p space shows 
that for 0.60 < p < 1.05 a.u. the line dimensions of the FS 
along [001] are much larger for p || TK than for p \\ TM, 
which is caused by the lack of electrons in the 3rd and 4th 
bands around the L point and the lack of holes around 
the H point. The anisotropy of the Fermi momentum, 
i.e. p™ > Pp K , leading to very small holes around the 
£ point, is also observed. All subtle features of the the- 
oretical anisotropy are reproduced by the experiment in 
detail. However, its magnitude is diminished, which indi- 
cates that e-e correlation has an anisotropic effect on the 
momentum density. With the example of the anisotropy 
Prjkip) ~ Ptm(p)i ^ i s shown how symmetry properties 
can enhance the accuracy of the reconstruction. 

After folding the reconstructed densities p 001 (p z ,Py) 
into the reduced zone scheme, the 2nd zone holes and 3rd 
zone electrons could be clearly recognized around the T 
and K points, respectively. Moreover, the results both 
in p and k spaces clearly show the absence of the FS 
elements around the L point in the 3rd and 4th bands. 
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